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This article presents a calculation of the mean electromotive force arising from general small-scale magneto- 
hydrodynamical turbulence, within the framework of the second-order correlation approximation. With the goal 
of improving understanding of the accretion disk dynamo, effects arising through small-scale magnetic fluctu¬ 
ations, velocity gradients, density and turbulence stratification, and rotation, are included. The primary result, 
which supplements numerical findings, is that an off-diagonal turbulent resistivity due to magnetic fluctuations 
can produce large-scale dynamo action - the magnetic analogue of the “shear-current” effect. In addition, con¬ 
sideration of a effects in the stratified regions of disks gives the puzzling result that there is no strong prediction 
for a sign of a, since the effects due to kinetic and magnetic fluctuations, as well as those due to shear and 
rotation, are each of opposing signs and tend to cancel each other. 


I. INTRODUCTION 

Explaining the amplification of magnetic fields with cor¬ 
relation lengths larger than the underlying fluid motions has 
proven to be a fascinating and rich problem in astrophysics. 
From the early days of mean-field dynamo theory it has been 
well known that the presence of fluid helicity enables such be¬ 
havior [1, 2]. This is the so-called a effect, where the small- 
scale turbulence creates an electromotive force (EMF) 

{uxb) 

that is proportional to a large-scale magnetic field, £ = aB, 
leading to exponential instability in the kinematic regime. 
While this simple a effect is now well established and reg¬ 
ularly observed in simulations, a variety of complications ex¬ 
ist in explaining observations. For one, in some situations - 
for instance, the inner regions of accretion disks - there is 
no reason to expect a helical flow and symmetry arguments 
demonstrate that a = 0, yet dynamo action is still observed 
in numerical experiments [3, 4]. Fess obviously, nonlinear 
effects caused by the fast build up of small-scale fields can 
“quench” a dynamos before significant mean-field amplitudes 
are reached [5, 6]. Since the effectiveness of this quench¬ 
ing increases with the Reynolds numbers, it remains unclear 
whether mean-field theory is able to explain the observed field 
amplitudes in the nearly dissipation-free plasmas prevalent in 
astrophysical environments. For these reasons, it is interesting 
to consider other possibilities for mean-field dynamo action, 
in particular the effects of velocity gradients and strong ho¬ 
mogenous magnetic fluctuations. 

In this paper, we present a very general theoretical examina¬ 
tion of different mean-field dynamo effects, within the second- 
order correlation approximation (SOCA). In particular, we in¬ 
clude the effects of specified large-scale velocity gradients, 
rotation, density and turbulence stratification, helicity, and a 
bath of strong small-scale magnetic fluctuations (treated in 
the same way as the velocity fluctuations). For our primary 
inspiration in this work - the accretion disk dynamo - each of 
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these effects can be important in some way, and this will also 
be the case in a wide variety of other astrophysical scenarios. 
Of particular note is the presence of homogenous magnetic 
fluctuations, which have not been included in most previous 
theoretical mean-field dynamo investigations (but see, for ex¬ 
ample, Refs. [7-10]). These should be generically present, 
at a similar level to velocity fluctuations, in magnetohydrody¬ 
namic (MHD) turbulence above moderate Reynolds numbers, 
due to small-scale dynamo action. While SOCA itself cannot 
capture the small-scale dynamo, by assuming the presence of 
the magnetic fluctuations we can compute expected changes 
to the EMF, in particular whether a small-scale magnetic field 
might suppress, or enhance, kinematic dynamo effects. 

The most important result presented here is an analytic 
confirmation of our numerical work related to the “mag¬ 
netic shear-current effect” [11, 12]. Generically, this type of 
dynamo is non-helical, driven by the interaction of an off- 
diagonal turbulent resistivity with a mean shear flow [13-16]. 
Some controversy has surrounded the kinematic version of 
this effect, since following early work [14—16], others found 
that the crucial transport coefficient rfy^ had the incorrect sign 
to promote dynamo action [17-20]. Here, we show that the 
magnetic version of this effect is much more robust and of the 
correct sign - not only is its magnitude substantially larger 
than the kinematic effect, but a variety of calculation meth¬ 
ods agree on this: SOCA, the spectral t approximation [8], 
quasi-linear theory [11, 20], and perturbative shearing wave 
calculations [21]. With this array of other calculations, we 
feel that SOCA calculations are important, not because they 
should be more accurate than other methods, but because they 
are simple, have a well-understood range of validity, and al¬ 
low exploration of expressions across a range of parameters 
(e.g., magnetic Prandtl number). This final consideration is 
notable since it provides the researcher with some indication 
of the robustness of a given effect, for instance by noting if 
the sign a given transport coefficient is particularly sensitive to 
slight changes in parameters. Finally, all of our results related 
to rjyx have been confirmed through direct numerical simula¬ 
tions [11, 12]. Most important is the measurement of a marked 
decrease in rjyx after saturation of the small-scale dynamo in 
sheared turbulence, accompanied by excitation of a coherent 
mean-field dynamo [12]. 

Turbulence and density stratification is invariably signifi- 
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cant in astrophysical scenarios, including in accretion disks 
away from the central plane of the disk. With this application 
in mind, we also apply our results to the case of stratified ro¬ 
tating turbulence with strong velocity shear, considering the 
resulting a effects. We find that for a Keplerian (or more gen¬ 
erally, anticyclonic) rotation profile, the contributions from 
shear and rotation, and those from kinetic and magnetic fluctu¬ 
ations, are each of opposite signs. The dominant contribution 
will depend strongly on the magnetic Prandtl number Pm, as 
well as the relative intensities of magnetic and kinetic turbu¬ 
lence. This is confusing in light of the beautifully coherent 
“butterfly diagrams” that are often seen in stratified accretion- 
disk simulations [22-25], which would suggest a robust neg¬ 
ative value for ayy. We note that the contributions to these a 
effects from velocity shear are at least as strong as those from 
rotation and should not generally be neglected. 

The structure of our calculation almost identically follows 
that of Radler and Stepanov [18] (hereafter RS06), with the 
additional effects of magnetic fluctuations, density stratifica¬ 
tion (within an anelastic approximation) and net helicity. The 
inclusion of such a variety of physical effects leads to a rather 
prodigious number of terms, and we have used the VEST 
package [26] in Mathematica to carry out the bulk of the cal¬ 
culations. We start, in Sec. II, by outlining the setup of the 
calculation, including the most general form of £ allowed by 
the symmetries of the problem, as well as the relation of the 
transport coefhcients in Cartesian domains with velocity shear 
to this general form. We also give the perturbation expansion 
used, which is a generalization of that in RS06 to include mag¬ 
netic turbulence at lowest order. In Sec. Ill, we outline the pro¬ 
cedure used in the calculation itself, skipping many details for 
the sake of brevity. Particular focus is placed on the unstrati¬ 
fied shear dynamo - especially the magnetic shear-current ef¬ 
fect - in Sec. IV, while the stratified a effect is examined in 
the same geometry in Sec. V. Readers interested primarily 
in the application of calculated coefficients to disk dynamos 
may wish to skip directly to these sections. Due to the length 
of algebraic expressions, the full set of transport coefficients 
is given in Appendix B. 


II. FUNDAMENTALS OF MEAN-FIELD 
ELECTRODYNAMICS 

Our starting point, common to most mean-field dynamo 
calculations, is the system of compressible MHD equations, 

§ + V ■ ipUr) - 0, 
ot 

dUr 

p - \- piJJj • V) + — Bj • V Bf 

ot 

-t- V ■ [pvCVUt + (VUTy) + p^dijV ■ Uj] -H <T„, 

^^Vx(UtxBt) + + a-b, 

ot 

( 1 ) 


Here Uj and Bj are the full velocity and magnetic fields, v is 
the kinematic viscosity, is the bulk viscosity (this will not 
contribute), and fj is the resistivity. We have included the ef¬ 
fects of rotation through a mean Coriolis force (2pQ xUj) 
in the momentum equation. Before calculating transport co¬ 
efficients from Eq. (1), we shall apply an anelastic approxi¬ 
mation [27, 28], assuming nearly incompressible fluctuations 
with V ■ ipu) - 0 [see Eq. (2)]. This allows low-order ef¬ 
fects due to a mean density gradient to be retained, while still 
preserving most of the simplicity of an incompressible calcu¬ 
lation. 

Mean-field dynamo theory [1, 2] involves splitting fields 
into a mean and fluctuating part; 

Ut^U + u, Br^B + b, ( 2 ) 

with U = (Ut), B = (Bj). The averaging operation (•) should 
filter out small scales and satisfy the Reynolds averaging rules 
(later in the manuscript we will specify (•) as a horizontal spa¬ 
tial average). Applying (■) to the induction equation leads to 
the well-known mean-field induction equation 

dtB x(UxB) + V x8 + vaB. (3) 

where £ = (m x ft) is the electromotive force (EME). The goal 
of mean-field theory is to calculate £ as a function of B and 
other parameters in the problem (i.e., U, Q, Vlnp and the 
small-scale turbulence statistics), thereby closing Eq. (3). If 
£ (B) is such that a small magnetic field will be reinforced by 
the small-scale turbulence, a dynamo instability results. 

Before commencing with a full calculation of £, it is worth 
examining the symmetries of the problem. Assuming scale 
separation between the mean and fluctuating fields, we can 
Taylor expand the EME as 

£/' — + . . . ( 4 ) 

where we use the Einstein summation convention and the 
comma denotes a derivative. The tensors aij and bijk are the 
transport coefficients determined by the turbulence. In keep¬ 
ing with the separation of scales assumption, we shall consider 
linear B fields (B),- = Bi + BijXj, velocity fields (17),• = UijXj 
and density p = po H- po ■ V Inp (the constant velocity part 
can be removed by Galilean transformation). As in RS06, 
to cleanly separate different dynamo effects into scalar coef¬ 
ficients, it is helpful to split VU and VB into symmetric and 
antisymmetric parts, 

Bij = (VB)« -Afj^ (VB)« - ^-e^jkJk, (5) 

where and are the symmetric and antisymmetric 

parts of Uij and By, W - V x 17 is the background vorticity and 
J - VxB is the mean current. Due to the assumption V-17 = 0 
in our calculation, we have implicitly assumed 17 • Vp = 0, a 
requirement that could easily be relaxed if desired. 

We consider general inhomogenous background turbulence 


V ■ l/r = 0, V ■ B, = 0. 
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in both u and b, modified by mean velocity gradients, rotation 
and density stratification. The density stratification is assumed 
to be aligned with the turbulence stratification in the direction 
g, but we allow their magnitudes and signs to differ; that is, 
defining 


Vlnp^Xpg^ Vlnw ^ XuS, Vlnl? ^Xhg, ( 6 ) 

(where we allow;t'p + xn + Xh- 

For completeness, we include both non-helical and helical 
contributions to the turbulence [29] but neglect the effects of 
inhomogeneity on the helical part [30]. We assume that the 
EMF due to the background turbulence vanishes, {u xb}^ - 
0. Such a B independent contribution could be important in 
some situations (see, for example, Yoshizawa and Yokoi [31]) 
and the method applied here can be used to calculate well- 
known effects of this type if desired, for instance the cross- 
helicity effect [32]. In addition, we do not calculate the com¬ 
ponents of the Reynolds stress, which would force a mean- 
field velocity U. This is not justified for any particular reason 
other than our primary interest in the magnetic field dynam¬ 
ics. While it is possible that there are important interactions 
between U and B that lead to other instabilities [33], we leave 
their systematic study to future work. 

A careful consideration of the symmetry properties of the 
system leads to the general representation of £ in terms of a 
set of scalar transport transport coefficients (see RS06 for a 
full explanation) 

£ = - xB- x B 

-af\g • Q)B - af[{g ■ B)n + (B • Q)#] 

-a^Pig ■ W)B - a''^\(g ■ B)W + (B ■ W)t] 

tr^ ^i,Silm^ljgm ^jlmb^ligm)Bj 

_(y(0) + y(n)^ X Q + y(W)^ ^ ^ y^^^Dijgj) X B 

- . (VB)« - iK^^^SijkDtr iVB)f (7) 

Here we have conformed to the sign conventions in RS06 and 
use the Einstein summation convention. The subscript -h de¬ 
notes a coefficient that is only allowed by the helical part of 
the turbulence, while all other coefficients arise only through 
the nonhelical part. In addition, since we assume small-scale 
fluctuations in both u and b, we further split each transport co¬ 
efficient into these contributions; e.g., 

Since we work with SOCA in the linear regime (where B is 
small), these are always additive and can be calculated sepa¬ 
rately from the u and b turbulent contributions. 


A. Cartesian domains 

In Sec. IV we shall give specific results for the numeri¬ 
cally convenient Cartesian shear dynamo with nonhelical, un¬ 
stratified background turbulence. This is essentially a gen¬ 
eralization of the unstratified shearing box that is often used 


in accretion-disk simulations. In this case, mean fields de¬ 
pend only on z, U - -S xy (giving W - -S z), = Qz and 

the mean-field average is defined as an average over x and y, 
(•) = (LxLyY^ f ■ dxdy. The mean-field equations simplify to 

diB^ = -rjy^dlBy + ijyydlB^, 

d,By ^ -SB,- J],ydlB, + T]„dlBy, (8) 


where the rjij are defined to be the relevant components of bijk 
that are nonzero for the chosen average and mean field. Eor 
Bi — BiQe‘^^e^‘ a coherent dynamo is possible if 


F = k\-ST]y, + k 


Tlxyilyx + -idxx - rtyy) 


1/2 


- k^iBxx + dyy) 


(9) 

has a real part greater than 0. One can neglect the term mul¬ 
tiplying k^ in the square root in Eq. (9) since S is presumed 
to be large compared to all transport coefficients. This gives 
t] 2 iS < 0 as a necessary condition for instability. Computing 
the relationship between Eq. (7) and Eq. ( 8 ) shows 



^(n)_ 


( 10 ) 


and T],, - rjyy - Note that Eq. (9) only describes the 
growth due to a coherent dynamo process and fluctuations in 
a or 77 that arise in any finite system can cause a dynamo in and 
of themselves [11, 34, 35]. We shall specialize to the Cartesian 
case in Sec. IV and keep U general for the calculation of the 
transport coefficients listed in Eq. (7). 


In Sec. V we give results specific to the case of stratified 
sheared rotating turbulence. This is motivated by consider¬ 
ation of the upper (or lower) portions of an accretion disk. 
Again, mean fields depend only on z, U = -S xy, Q - Qz, 
and g - z. We neglect off-diagonal resistivity contributions 
and use rj,, — Tjyy — 0'^^ ■ The mean-field equations simplify 
to 


d,B, = -ay,d,B, - ayyd.By + BIb, 

d,By = -SB, + a,yd,B, + a„d,By +0^^dlBy. (11) 

With a,y - -Oy,, considering B, = one obtains the 

growth rate 

T -{ikSayyl2 + 0ayya,^^ +ika,y-00^\ ( 12 ) 

Again, S is presumed large in comparison to all transport co¬ 
efficients, so we see that any nonzero Oyy can lead to insta¬ 
bility at sufficiently long wavelength. Of course, in practice 
there will be a minimum k possible in the system, particu¬ 
larly since ayy arises from a stratification, so a finite Oyy will 
be necessary to overcome the turbulent resistivity. The co¬ 
efficients in Eq. (11) are related to those in Eq. (7) through 
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axy - -Oyx - and 

Qyy - S - Qaf^\ 

Oxx = S - Q af\ (13) 


B. Perturbation expansion to describe tbe fluctuations 


For the calculation of £ we use the second-order corre¬ 
lation approximation (SOCA), which involves solving linear 
equations for the fluctuations by neglecting third-order and 
higher correlations. As such, this is rigorously valid only 
at low Reynolds numbers where dissipation dominates over 
nonlinearities for the fluctuations (SOCA can also be valid in 
the small Strouhal number limit [Eq. (31)], see Brandenburg 
and Subramanian [36] for a more thorough discussion). In 
addition, we choose to include the shear, rotation and den¬ 
sity stratification perturbatively [16, 18], considering only the 
linear response of transport coefhcients to these effects. An 
analytic calculation with shear included at zeroth order can 
be found in [20], and some examples of calculations that in¬ 
clude nonlinear contributions from other effects can be found 
in Refs. [8, 9, 27, 28, 37]. In a very general calculation, Pipin 
[10] nonlinearly includes all effects discussed here (although 
the approach, the “minimal t approximation,” has a some¬ 
what unknown range of validity). We have also computed the 
magnetic dynamo transport coefficients with non-perturbative 
shear and rotation using statistical simulation in the shearing 
box [11]. 

Following Rudiger and Kichatinov [28], Kichatinov and 
Rudiger [27], and Rudiger [37], we start by making an anelas- 
tic approximation to the full compressible equations, V-(pu) - 
0. This should be valid for weakly compressible turbulence 
and allows the inclusion of a weak density stratification into 
the problem, which is important in a wide variety of mean- 
field dynamos. We shall assume that the large-scale flow is in¬ 
compressible, since our primary application is to shear flows. 
It is then more convenient to work in terms of the small-scale 
momentum [27, 28], m = pu, since the calculation for m pro¬ 
ceeds in a similar manner to the incompressible case. 

In retaining both strong homogenous velocity and mag¬ 
netic fluctuations, denoted uq (or mo) and bo respectively, 
we must treat the momentum and induction equations on the 
same theoretical footing. We start from Eq. (1) by split¬ 
ting into mean-field and fluctuation equations, applying the 
anelastic approximation followed by the change of variables 
Uq - ntolp. We then linearize the small-scale equations and 
expand m = mq + b - bo + ft® -H ft*'^..., to 

perturbatively find the change to the background turbulence 
caused by the shear, rotation and stratification. This leads the 
SOCA equations that will be used to calculate all transport 


coefhcients: 

(dt - vA) m® = - (two ■ VI/ -H 1/ ■ ViMo - (gp ■ l/)»io) 

-Vp® - 2n X mo H- (fto ■ VB + B ■ Vft®) - vgp ■ V/mq 
(d, - VA) -VU+U- - {gp ■ 

- 2n X m® H- (ft® VB + B Vft®) - vgp ■ Vm® 
{dt - i]A) ft® = [(gp • mo)B -mo -VB + B ■ Vmo 

- {gp ■ B)mo\ + bo-VU -U ■ Vfto, 

{d, - rjA) ft<') = p-' [(gp • m®)B - m® VB + B Vm® 

- (gp ■ B)m®] H- ft® -VU-U- Vft®, (14) 

along with divergence constraints for each m®, ft®, and 
ft'^^^ Here gp s ;^fpg and we have neglected second derivatives 
of U and p, as well as products of VB with Xp [these contri¬ 
butions should vanish in the transport coefhcients, since the 
Eq. (7) illustrates that there is no contribution to the resistiv¬ 
ity due to g at linear order]. In addition, we shall neglect any 
terms that involve quadratic products of U, £1, and Xp (c-g-^ 
(gp ■ V)mQ), and expand all terms to linear order to take the 
Fourier transport of Eq. (14) (see App. A). 

While it may seem surprising that one requires terms two 
orders higher than mo and bo, it is straightforward to see that 
only considering and b^^'’ will not lead to contributions 
to £ that depend on products of B with U or £l (these are 
the interesting terms in the dynamo, describing the effect of 
rotation or velocity). With this in mind, the EMF is calculated 
as 

&ij ^[uib^ = [p^^moiboj) + (p"^mo,ft®) -H (p^^maft®) 

+ + (p^^m®fto;) H- (p“'m®ft®) . (15) 

Despite the fact that all the terms in Eq. (15) give some con¬ 
tribution, there are also a large number of terms that con¬ 
tain quadratic products of Up, Q,, Xp^ or B, which are ne¬ 
glected. As is evident, with background turbulence in both u 
and ft there will be contributions to £ from the Maxwell stress 
{B ■ Vb + b ■ VB) that one would expect to be of a similar 
magnitude to the standard kinematic dynamo arising from the 
Lorentz force [V x (m x B)]. This choice of perturbation ex¬ 
pansion is the natural generalization of RS06 to the case with 
bo fluctuations (although note that in RS06 has become 
in our notation such that u and ft are treated on equal foot¬ 
ings). Our results for the kinematic dynamo {bo - 0) without 
density stratification agree with RS06 aside from a single nu¬ 
merical coefhcient (see App. B). 


III. OUTLINE OF THE CALCULATION OF £ 

Our calculation follows the methods and notation in RS06 
and a full explanation is given there. Here we give a very brief 
outline, in particular the choices involved, with final results 
given in Appendix B. We have carried out the entire calcula- 
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tion in Mathematica using the VEST package [26] to handle 
abstract tensor manipulations using the Einstein summation 
convention. 

The two-point correlation of two fields v and w is defined 
as 


{x\tVX 2 ,t 2 ) = {vi{x\t\)wj{X 2 ,t 2 )) ■ (16) 

It is convenient to write such quantities in the variables 

= (jti -H X 2 ) /2, r - x\- X 2 , 
r = (fi H-f2)/2, t - t\ - t2, (17) 


giving 

( 18 ) 

One then Fourier transforms in the small-scale variable r to 
obtain 


{R, T\r,t)^ j dk doj {R, T\ k, co) ( 19 ) 

with 

^l'"HR,T-,k,u) = J dKdil ([v]+ (20) 

where v - v{k, a>) and w - w(k, co) denote the Fourier trans¬ 
forms of V and w, and we use the [•] j, notation of RS06, 

[f(k,co)]^^f(±k + K/2,±aj + D./2). (21) 

As in RS06 we shall calculate 


3ij (R, T; 0,0) = j dkdto Sij {R,T-,k,u)) 

= J dKdQdkdco{[p-^mi]+[bi]-)e‘^ ‘^^‘^'^ 

= f dKdQdkdupo^{[mi-igpjdkmiUbi]-)e‘‘^‘^^‘^^, 

( 22 ) 

setting /?, r —> 0 only after extracting the coefficients of B, 
and Bij (i.e., the transport coefficients ay, bijk). 

With these notations defined, the starting point of the cal¬ 
culation is the substitution of the linear forms for U, p and B 
and into Eq. (14), followed by a Fourier transform. This leads 
to Eqs. (A2)-(A5). One then substitutes mf^ and Vf' into 
and yf’ to form explicit expressions for m, and hi in terms of 
Mo; and hiQ. Defining 

niij = ([mo,]+[mo2]-), 
bij^{[boi\Ahj]-), 

to specify the statistics of mq and b^, this allows one to form 
Eq. (15) in terms of my and bij, neglecting all terms that con¬ 
tain UijUrs, Uij^j, Q.iQ.r, UijXp, kliXp, (Vlnp)2, or any prod¬ 


ucts of Bi and By. Recall that we have assumed (mq^o) = 0, 
implying that all terms in the expansion of £y contain B, or 
By. In keeping with the expansion to linear order in back¬ 
ground quantities, it is necessary to expand [f(k)]± to first 
order in K in those terms that contain Bi (i.e., a coefficients). 
These lead to terms involving the gradient of the turbulence 
intensity. Note that [f(k)]± —> f(±k) for resistive terms (coef¬ 
ficients of By). Some useful identities in the above procedure 
are given in RS06 Eqs. (33)-(35), which are needed to remove 
d/dki derivatives from mq, and boi. Similarly, we apply the 
identities 


kiMij = kiiTiji = —mji (23) 

(and similarly for hjf), which arise from the divergence con¬ 
straints on m, and hi. 


Extracting the coefficients of B, and By in the expression 
for £, = Sijk&jk (0,0), at this stage we have large integral ex¬ 
pressions for oy and foy*. in terms of my and foy and their spa¬ 
tial derivatives [for example, RS06 Eqs. (39)-(40)]. Without 
further interpretation, such expressions are nearly useless, and 
it is helpful to insert explicit forms for my and hij. Assuming 
isotropy in the limit of vanishing mean flow and rotation, we 
insert 


bii = - 


- isiji^Hu (k, o )), 
-^2^kiKj-kjK,) 


Wm (K; k, to) 


Wb (K-,k, to) 


- iSiji-j^Hb(k,(o), 


(24) 


where k - |k,|. Here Wm,h represents a non-helical part and 
Hm.h a helical part of the background turbulence [27, 28]. This 
form for Wm is particularly convenient since it can be shown 
that to first order in the scale of density variation 

Wm (x', k, to) - p^(x)Wu (Jt; k, to ), (25) 


where Wu (x; k, co) is a similar function specifying the statis¬ 
tics of u and Wm(x',k,u) - f dKe‘^'^Wm(K;k,a>) [27]. In 
this way. 


VWm(x; k, oj) = g(2xp + 2xu)Wm (x; k, to ), (26) 

separating the effects due to density and turbulence stratifica¬ 
tion. Similarly, for the magnetic fluctuations 

VWb(x; k, to) - 2gXbWb (x; k, to). (27) 


It transpires that all terms now depend on k only through k. 
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and all of the integrals can be substantially simplified using 


/ 


dkkikjfik) - -Sij I dkk^f(k). 


M UMJ m = * «,.) / dk m. 

(28) 


where the integrals over k on the right-hand side of Eq. (28) 
are taken from k = 0 — > oo. One then splits Uij and 
Bij using Eq. (5), putting £, in the form given by Eq. (7). 
One can straightforwardly read off the transport coefficients 



,..., as integrals of the form 


J dk doj 00^ (k, u) Hu,b (k, (d) 

= 4-71 

\ 'u,b 

j dk doj 00'^ (k, to) Wu,b (k, to) 


dk dto 00'^ (k, to) Wu,h (k, to). 


The full list of coefficients ... is given 

in App. B. 

Einally, it is possible to carry out the integrals of the form 
in Eq. (29) for a specific form of W and H, leading to explicit 
expressions for the transport coefficients in terms of the physi¬ 
cal parameters. A convenient form for examining expressions 
and plotting is the Gaussian W used in RS06, 


2AItc 

3 1 + 


(30) 


with a similar definition of Wb- With this choice, all integrals 
can be carried out explicitly without further approximation. 
As in RS06, we shall write such expressions in terms of the 
non-dimensional variables (and po) 

e = fimis/Mims, P = d^/vTc, q = A^lpTc, Pm = vlq. 

Re — Mrmsdc/v, Rm — Mj-ms/ic/r/, St — Mj-jnsTc/dc. (31) 


Here Pm, Re, Rm, and St are respectively the magnetic Prandtl 
number, the fluid Reynolds number, the magnetic Reynolds 
number and the Strouhal number, p and q are the ratio of dif¬ 
fusion times, d^/v and d^//;, to the correlation time Tc- Thus 
^ > 0 denotes the low conductivity limit, while q ^ oo de¬ 

notes a high conductivity limit (with a similar result for p 
and fluid diffusivity). A sufficient condition for the validity 
of SOCA (i.e., neglect of nonlinear terms in the correlation 
equations) is Rm <K 1 in the limit q —> 0, and St «: 1 in 
the limit q ^ oo, see Brandenburg and Subramanian [36] and 
Radler and Stepanov [18] for more discussion of these validity 
regimes. In addition, e require Uij and Q, be a small pertur¬ 
bation to the background turbulence. In practice, we shall use 
these non-dimensional variables [Eq. (31)] for plotting trans¬ 
port coefficients. 

We have carried out the full sequence of steps detailed 
above in Mathematica using the VEST package [26] to en¬ 


able straightforward manipulation of tensors in index nota¬ 
tion. This has the obvious advantage of handling the very long 
expressions with ease and making the calculation straightfor¬ 
ward to generalize or modify. The sequence of steps is es¬ 
sentially the same as that detailed above. We first define m^'^\ 
^(0), ^(1)^ insert and ft® into and ftP\ then 

only later remove products that are quadratic in t/y, Q, orXp- 
It is then straightforward to define [■]j. operators, their associ¬ 
ated product rules, and methods to in expand in K. This allows 
the construction of the entirety of £ in one step. Insertion of 
the explicit forms for vy and fty [Eq. (24)] and the partial inte¬ 
gration using isotropy [Eq. (28)] is easily carried out using re¬ 
placement rules. Einally, we decompose products of By with 
Uij, il and g into the form given in Eq. (7), allowing the co¬ 
efficients listed in App. (B) to be straightforwardly extracted 
from the total expression. Einally, if so desired, these can be 
directly integrated with the specific form of W [Eq. (30)] by 
carefully substituting the dimensionless variables [Eq. (31)] 
and using Mathematica"s native Integrate function. For the 
interested reader, we include the full calculation notebook in 
supplementary material. 


A. Agreement with previous works 

Our results agree with related works of other authors in 
special limits, including those utilizing different calculation 
methods. As discussed throughout the work, all results of 
RS06 are recovered in the limit V Inp = 0 [aside from one dis¬ 
crepancy, in This agrees with Rudiger and Kitchati- 

nov [17], many results of Pipin [10], including his magnetic 
contributions (see his App. B), as well as the quasi-linear 
methods in Sridhar and Subramanian [38] and Singh and Srid- 
har [20]. As is well known, there is a discrepancy between 
these kinematic quasi-linear results and those obtained using 
the T approximation [8, 16], possibly due to a change in sign of 
Pyx with Rm [19]. As seen in Eq. (32) of Pipin [10], his con¬ 
clusions regarding the kinetic and magnetic contributions to 
the shear-current effect (with rotation) are are similar to ours. 
Our results also compare favorably to previous works without 
velocity gradients, but including magnetic fluctuations. As 
expected, the helical magnetic a effect has the opposite sign 
to the kinematic effect, and there is no change to 0-^^ due to 
the addition of magnetic fluctuations. In addition, the signs 
of and (5^^^ agree with the t approximation calculation of 
Radler et al. [9] < 0, > 0, although there is not an 

exact cancellation at fi = ft as in Radler et al. [9]). 

The a effects arising through stratification and inhomo¬ 
geneity also show broad agreement with previous works. Be¬ 
cause of the linearity of the expansion in V Inp, U and Q., the 
density stratification contributes very little to the coefficients, 
aside from directly through VWm [Eq. (26)]. This means 
generally appears together with the turbulent gradient;^fu. The 
one exception to this is the “turbulent diamagnetism” term, 
which interestingly depends only on the turbulence gra¬ 
dient, not the density gradient, due to a cancellation (this is 
in agreement with Kichatinov and Rudiger [27]). Again our 
results without mean velocity broadly agree with the t approx- 
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imation magnetic turbulence results given in Radler et al. [9]; 
for instance, the fact that = - (y°) and the opposing 

signs of the rotational kinematic and magnetic diagonal a ef¬ 
fects {a^)u,b, with |(q'^)„| > |(Q'^)i| (although we see a strong 
dependence of these parameters on Pm; see Sec. V). 


IV. SPECIFIC RESULTS FOR UNSTRATIFIED SHEAR 
DYNAMOS 


In this section we discuss the results pertinent to our pri¬ 
mary motivation for this work, the shear dynamo in a Carte¬ 
sian box. As shown in Eq. (8), in this geometry with a hor¬ 
izontal mean-field average, the number of transport coeffi¬ 
cients reduces significantly. We are particularly interested in 
the sign of the rjyx coefficient, which should be most impor¬ 
tant for dynamo growth due to its coupling with the shear 
[Eq. (9)]. Here we outline the contribution to Tjyx from ve¬ 
locity and magnetic fluctuations in the presence of shear, both 
with and without rotation. This geometry is particularly rel¬ 
evant for the central regions of accretion disks, where there 
is strong flow shear, stratification may be subdominant, and 
there is no obvious source of helicity in either velocity or mag¬ 
netic fluctuations [4]. 

Utilizing Eq. (10) and the results in listed in App. B, one 
obtains after some impressive cancellations 


(^yx)u 


/ 


32nk^Wu{k, ai) u/ff 

dcodk -;-, 

15 {if- + o)^) (v^ -H aP-) 


(32) 


/ 


idyx)b - I dojdk^nk^p ^Wb{k,oS) 


4^4 


15 (v^ + aP)^ 
2fi^ + fj^\P + 2u^fj^ + 

15 ( 77 ^ + aP) (v^ -H aPY' 

4tdPfjv 

15 (if- + aP)^ (v^ -I- aP), 


(33) 


i%x) 




,2^2 


dcodk 


64nk‘^Wu (k, oS) 

15 {ff + (jp)^ (v^ -H cP) 


(34) 



FIG. 1. Transport coefficients (%x)f (solid, blue), (rjyx)^ (dashed, 
blue), (Tiyx)l (solid, orange) and (r/yx)^ (dashed, orange) as a function 
of q for (a) Pm = 1, (b) Pm = 10, and (c) Pm = 1/10. Each coeffi¬ 
cient has been calculated using the form given in Eq. (30) for IV, and 
normalized by with the magnetic diffusion time, A^lq, held 

constant (equivalently = 1/^). (Note that this choice is necessary 
because the coefficients have different units, and is chosen purely for 
plotting purposes, since it reduces the variation of coefficients with 
q.) 


(dyx)h ^ - dai dk 


Snk^p ' Wb {k, to) — l2aPP + 3aj^) 


15 (v2 + to^f 


(35) 


Here v = vP, rj - qP, integration over to is from -00 to 00 
and over k is from 0 to 00 . We have defined each coefficient 
such that 


Tlyx = s [{riyx)l + {qyx)l] + ^ [{riyx)^ + (7rx)“]. (36) 

to keep all signs consistent. Recall from Eq. (9) that with our 
definition of S , qyyS < 0 is required for a growing dynamo 
(note that this is the reverse of RS06). Eor Keplerian rota¬ 
tion, Q. -2S13, since vorticity and rotation are opposite (i.e.. 


anticyclonic) when S and Q have the same sign. 

Let us first examine the coefficients for a kinematic dy¬ 
namo, i.e., with strong homogenous velocity fluctuations [the 
coefficients (qyx)u, Eqs. (32) and (34)]. Eirstly, we note that 
the contributions from S and Q. have identical forms, and that 
the integrands are positive definite [39], see Eig. 1. Thus, as is 
well known, we see that iqyx)l , the “shear-current effect,” has 
the incorrect sign for dynamo action within this quasi-linear 
approximation. Although the basic fix J effect (also known 
as the Radler effect) is well known, the explicit calculation 
of transport coefficients including shear and rotation seems 
to have been mostly ignored, although there is much discus¬ 
sion in early literature on the subject (e.g., Krause and Radler 
[2], Moffatt and Proctor [13]). Given the identical forms of 























Eqs. (32) and (34), we can immediately write down the result 

(/7.yx)« = (5-2Q)H, (37) 

where H is the (positive) integral in Eq. (32). Thus, we find 
that the addition of Keplerian rotation (Q = 2S /3) (as relevant 
to turbulence in accretion disks for example), will change the 
sign of rjyx to slightly negative and a coherent dynamo insta¬ 
bility should be possible. Indeed, this is seen in our recent 
simulation work [ 11 ], where we observe increasing coherency 
and a larger growth rate as the rotation is increased in the an- 
ticyclonic direction. 

Turning to the coefficients for magnetic fluctuations we find 
the interesting possibility of a magnetically driven dynamo. 
In particular, as shown in App. (C) and Eig. 1, the coefficient 
ij]yx)h is consistently negative and generally larger than the 
other contributions. This implies that a dynamo can be excited 
by magnetic fluctuations, themselves presumably arising from 
a small-scale dynamo process, or perhaps an MHD instability 
of some sort. Since the small-scale dynamo is usually consid¬ 
ered harmful to mean fields [5], this is an interesting possibil¬ 
ity - a build of magnetic noise on small scales may cause a 
coherent large-scale dynamo to develop. The addition of rota¬ 
tion renders the effect of magnetic fluctuations more complex, 
and no simple result seems possible. In particular, the sign of 
the (riyx)^ coefficient depends on the parameters, and is gen¬ 
erally negative for large v, rj and positive at lower dissipation, 
although smaller in magnitude than {riyx)\ This change in sign 
is also seen in quasi-linear calculations [ 11 ]; however, given 
that the quasi-linear approximation becomes less valid in this 
limit, it would be unwise to draw any conclusions about the 
high-Rm limit from this behavior. 

Einally, we note the possible relevance of this dynamo to 
the central regions of accretion disks. In self-sustaining tur¬ 
bulence simulations in this geometry, magnetic fluctuations 
are generally substantially stronger than velocity fluctuations. 
Such conditions seem ideal for excitation of a coherent dy¬ 
namo driven by the magnetic shear-current effect. We note 
that cyclic behavior, as often observed in self-sustaining sim¬ 
ulations [4, 25], seems to be quite generic in the nonlinear de¬ 
velopment of the magnetic shear current effect, and we have 
observed this in low-Rm simulations with a forced induction 
equation [11]. In addition, it is worth noting that Lesur and 
Ogilvie [4] concluded that rjyx was the primary dynamo driver 
from analysis of their numerical simulations. While more 
work is obviously needed to explore this possibility in detail, it 
seems reasonable to conclude that the magnetic shear-current 
effect is playing a fundamental role. 


V. SPECIFIC RESULTS FOR STRATIFIED ACCRETION 
DISKS 


Self-sustaining turbulence simulations in this geometry (for 
instance with shear-periodic boundary conditions in the radial 
direction) exhibit a very coherent dynamo, with quasi-time- 
periodic behavior in By and creating a “butterfly diagram” 
[23, 24]. Large-scale magnetic structures are seen to emanate 
from the central portion of the disk, migrating upwards into 
the lower density regions and becoming more intense as they 
do so [25]. This migration behavior would be characteristic 
of a dynamo driven by Uyy above and below the mid-plane: 
as shown in Eq. (12), growth of this type of “aw” dynamo is 
always accompanied by dynamo waves since E is complex. 
Note that a negative imaginary part of E is required for up¬ 
wards migration of mean-field structures with g - z. This 
occurs for Uyy < 0, a^y < 0, {Uyx > 0) [40]. 

Utilizing Eq. (13) with the results listed in App. B, and set¬ 
ting Pm = 1 here for simplicity, one obtains. 


r Wu(k, u)\^ (5\^ + cjA 

(ayy)l = ^nxp-u / dojdk -—--, (38) 

J 15 (v2 H- OJ^y 


(^yy)b - / dcodkp ^Wb(k,u)k' 


2 7v‘*-4wV-3w4 


15 (v^ H- w^) 


(39) 


r k^Wuik, u)v^ (v^ + 5w^) 

(«w)“ = -64;ryps / dojdk -—- -3 -, (40) 

J 15 (v^ + ax) 


„ r Wbik, u)kW (a? - 3v2) 

(ayy)b = -647r^5 / dujdk --- ^-3 -(41) 

J 15 (v 2 -H u^y 

Einally, for the off-diagonal component, 7 *°^ = Uxy = -Oyx, 
one has 

/ (Oh ^ f ^ k^Wu{k,u)fj 

(7 X = 47r^s / doidk—— -—, (42) 

J 3 ( 77 ^ -H 

Here we use the notation= |Vln(pM)|, and again signs are 
defined such that 

Oyy — S |^(Qy)i)„ + Xyy)b ] + ^ (^yy)b ] ’ (44) 

for anticyclonic rotation, e.g., Keplerian rotation is O = 2/35. 

It is first worth noting the sign of each coefficient given in 
Eqs. (38)-(43). W\i\\ Xpu, AS > 0 it can be shown easily from 
the above expressions that 


In this section we briefly outline how our results apply to 
stratified sheared rotating turbulence. Our primary motivation 
is consideration of the upper and lower regions of accretion 
disks, where the turbulence is stratified in density and inten¬ 
sity by the vertical gravity, perpendicular to the velocity shear. 


{ayy)l > 0, (ayy)l < 0, {ayy)^ < 0, {Uyy)^ > 0. (45) 

(Note that for the b components, it is necessary to integrate by 
parts over w, see App. C). The relations in Eq. (45) appear to 
also hold for Pm 4 ^ 1 (although we have a proof of this only 
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FIG. 2. Transport coefficients (ayy)^ (solid, blue), (dashed, 

blue), (a„,)j (solid, orange) and (ayy)l (dashed, orange) as a function 
of q for (a) Pm = 1, (b) Pm = 10, and (c) Pm = 1/10. Each coeffi¬ 
cient has been calculated using the form given in Eq. 30 for IV, and 
normalized hy with the magnetic diffusion time, A^jq, held 

constant (equivalently = llq). The dotted (black) curve in each 
plot shows the total a,,,, with equal kinetic and magnetic turbulence 
levels for Keplerian rotation, Q. = 2/35 [Eq. (44)], to illustrate the 
variability in these predictions. 


for the Q coefficients). This consistent difference in sign be¬ 
tween contributions is rather inconvenient for the application 
of SOCA results to stratihed accretion disks. Since one ex¬ 
pects <0,Xb<^ (although possibly > 0) [24, 41], we 
are left with the situation where not only do the a effects due 
to M and b partially cancel, but also those due to rotation and 
velocity shear! What’s more, as shown in Fig. 2, the relative 
contribution of each depends strongly on Pm. In particular, 
we see a dominance of {ayy)u over {ayy)b for Pm > 1, but this 
can reverse at low Pm. Similarly, the relative contributions 
due to velocity shear and rotation for the magnetic effect vary 
substantially with Pm, although the effect of shear seems gen¬ 
erally more substantial. While the ratio of kinematic shear 
and rotation contributions may be somewhat more robust, the 


two are roughly equal in magnitude, {UyyY^ - 

will approximately cancel for Keplerian rotation. Finally, it is 
worth noting that to complement these uncertainties, the signs 
of 7 *^°^ seem to predict the opposite held migration pattern to 
the upwards transport seen in simulation. In particular, for 
Xb < 0, Xit > 0, the kinematic and magnetic contributions 
both enforce > 0, leading to Im F > 0. However, in our 
use of the anelastic approximation, buoyancy effects are not 
included and these would be expected to change this aspect of 
the calculation substantially [40, 42, 43], potentially through 
large-scale instability [44]. 

Where does this leave us for understanding the dynamo 
in stratihed accretion disks? We see that aside from per¬ 
haps the transport term claims that SOCA predictions 
are incorrect for the stratihed regions of accretion disks are 
unfounded. More accurately, one could say that SOCA pre¬ 
dictions themselves are completely inconclusive, even in the 
kinematic regime, since each contribution - kinematic, mag¬ 
netic, rotation, and velocity shear - has a tendency to cancel 
its partner. Such uncertainty seems at odds with the robust dy¬ 
namo “butterhy diagram” seen across a wide variety accretion 
disk of simulations. 

Of course, one possibility is that the SOCA calculation car¬ 
ried out here, keeping only the linear contributions due to O, S 
and stratihcation, is not up to the task of calculating these co¬ 
efficients, and in reality there is a robust a effect. For instance, 
in Rudiger and Pipin [40], the authors hnd that ayy has the cor¬ 
rect sign (ayy < 0) for magnetic Huctuations in a compress¬ 
ible turbulence model for Keplerian shear and moderate Pm 
(this is the sign opposite to Eq. (39) but since their effect van¬ 
ishes in the incompressible limit, one should have no reason 
to expect agreement). Similarly, the calculations presented in 
Donnelly [45] go well beyond the accuracy of SOCA for the 
specific case of Keplerian shear through non-perturbative in¬ 
clusion of several extra physical effects; however, it is unclear 
from their (rather complicated) expressions whether the the¬ 
ory predicts a specihc sign for Oyy. While certainly feasible, 
it would seem a little bizarre that a behavior that appears so 
robustly in simulation could show so much variability across 
different calculation methods or rely on nonlinear behavior of 
transport coefficients with Q, 5, or the stratification. A va¬ 
riety of other possibilities might be imaginable, for instance 
a dynamo driven primarily by the magnetic shear-current ef¬ 
fect up to relatively far from the mid-plane (Sec. IV), with 
upwards transport above this caused by large-scale buoyant 
instability (not included here due to the anelastic approxima¬ 
tion). Another possibility could be that upwards held transport 
is caused by a small-scale magnetic helicity Hux [46, 47] from 
the central shear-current dynamo, which would create a (heli¬ 
cal) magnetic a effect. Such an process could look rather sim¬ 
ilar to a more standard a effect, although the basic cause of the 
dynamo would be entirely different [24]. Note that magnetic 
helicity Huxes have been found to be playing a signihcant role 
in unstratihed global MRI turbulence [48], providing some 
indication that such a process could be important. It is also 
worth noting that spatial variation in transport coefficients and 
quenching can lead to some interesting possibilities for dy¬ 
namo action [49, 50], and similar effects may prove important 
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at the boundary between the stratified and unstratified regions 
of disks. Overall however, it seems that the underlying cause 
for the “butterfly diagram” in stratified disks remains unclear 
and more work will be needed to arrive at robust mean-field 
models of the process. 

VI. DISCUSSION AND CONCLUSIONS 

In this work we have theoretically studied the dynamo in 
systems with mean velocity gradients, rotation, net helicity, 
and stratification, using perturbative calculations within the 
second-order correlation approximation. In addition to the 
standard kinematic dynamo, we have considered the possi¬ 
bility of a dynamo driven by small-scale magnetic fluctua¬ 
tions, as might arise from the small-scale dynamo or an in¬ 
stability. Our main finding is that an off-diagonal resistivity 
coupled to the shear can cause a dynamo instability in the 
presence of magnetic fluctuations. This effect - the mag¬ 
netic analogue of the “shear-current effect” [8, 16] - raises 
the interesting possibility of the small-scale dynamo enhanc¬ 
ing the growth of a large-scale field. In some sense, this pos¬ 
sibility is the reverse of large-scale quenching [5, 51]; rather 
than the small-scale magnetic fluctuations inhibiting the large- 
scale field growth, they could actively aid field generation, 
with large-scale growth eventually halting due to nonlinear 
changes to the transport coefficients, possibly influenced by 
secondary quenching effects [52]. 

Importantly, our prediction that the magnetic shear-current 
effect is able to excite a dynamo agrees with other transport 
coefficient calculation methods and simulations. In particu¬ 
lar, the T approximation predicts the linear magnetic effect to 
be much stronger than the kinematic effect (see Fig. 3 of Ro- 
gachevskii and Kleeorin [8]), just as was found in this work 
using SOCA (Fig. 1). In addition, agreement is found with 
quasi-Iinear calculations [11] (the magnetic version of the cal¬ 
culations in Singh and Sridhar [20]), as well as perturbative 
inhomogenous shearing wave calculations [21]. This suggests 
that the effect may be more robust than the kinematic shear- 
current effect and/or have less dependence on Reynolds num¬ 


bers. 

The work presented in this manuscript was primarily mo¬ 
tivated by gaining improved understanding of the fundamen¬ 
tal dynamo mechanisms in accretion disks. Consistent with 
the idea that two dynamo mechanisms might operate in disks 
[53], their inner regions seem well suited to be explained 
by the magnetic shear current effect [4] - magnetic fluctu¬ 
ations are generally stronger than kinetic fluctuations, rota¬ 
tion has the correct sign to enhance the kinematic dynamo, 
and the turbulence is essentially unstratified and nonhelical. 
Concurrent nonlinear direct numerical simulations of unstrat¬ 
ified shear dynamos in Cartesian boxes [11, 12] have con¬ 
firmed all results discussed in Sec. IV for the low-Rm regime 
[11, 54]. Firstly, we see a qualitative change in the kinematic 
dynamo with the addition of rotation, due to the change in 
sign of the rjyx transport coefficient [11]. Secondly, we observe 
the magnetically driven shear-current effect, both through di¬ 
rect driving of the induction equation [11], and at higher 
magnetic Reynolds number where magnetic fluctuations arise 
self-consistently though excitation of a small-scale dynamo 
[12]. The nonlinear saturation of these magnetically driven 
large-scale dynamos exhibits a pleasing resemblance to self- 
sustaining unstratified accretion disk turbulence simulations, 
with quasi-cyclic behavior of the large-scale By field. 

Less clear have been our findings regarding the a effect, 
as relevant to the stratified regions of accretion disks. In par¬ 
ticular, we find that a coefficients arising from rotation and 
shear, and those arising from kinetic and magnetic fluctua¬ 
tions, are each of opposite signs for anticyclonic rotation (Q 
and VxI7 antiparallel), and thus would tend to cancel. Further¬ 
more, predictions about which of these terms dominate (thus 
determining the sign of the total a effect), depend strongly on 
the magnetic Prandtl number and the relative levels of kinetic 
and magnetic turbulence. We thus conclude that perturbative 
SOCA calculations give no useful predictions regarding the 
primary driver of the so-called “butterfly diagram” pattern of 
large-scale field evolution seen in self-sustaining stratified ac¬ 
cretion disk simulations. Whether this is simply due to the 
inaccuracies of SOCA or there is some other more exotic ef¬ 
fect operating (e.g., a helicity flux [48]), remains to be seen. 
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Appendix A: Equations for b^°\ 6*'* in Fourier space 


Here we give the set of perturbation equations for u and b in 
Fourier space, as result from the Fourier transform of Eq. (14) 
The method is outlined in RS06, so we give very little detail 
here. Since we assume (/,■ (jt) = UijXj, p - po + XpSi^o and 
Bi (jt) = Bi + BijXj the Fourier transforms can be carried out 
exactly using Xkdibj - —dikbj — kidk^bj (where ' denotes the 
Fourier transform). We have also neglected products of Xp 
with Bij. In the momentum equations, the projection operator 
5ij - kikjlk^ is applied so as to remove the pressure. 

Defining, as in RS06, 




1 

ia> — vk^ ’ 


1 

ia> - rjk^ ’ 


the Fourier space equations are as follows, 


(Al) 




dniQi kikj 

-Uiinioi + - ivkrXpgrmm 


+ IV 


kikjkf 


. krD.r 


iy—j^Xpgrmoj + '^—j^SijkniQjkk + I 


ikyByboi 


- ikrB, 


kikj 


boi + Biiboi - Bikki 


dboi kikj 




boiB ji 


(A2) 




JJ (0) , TT I i . (0)jT ■ , - (0) 

-Uiim] + kfikki-^ + 2 — Uji- ivkrXpgrtn] 


kikjkf 


, • • J ' '' (0) , ,-,kr^r (0), , n jJO) 

+ iv—j^Xpgrnij + 2—j^Eijkm\ 'kk + ikyByb] 


kik j 


dbf^ kikj 


ib R ' I D iS^) _ Ji hr ^ _ O ^ ^ iu(Q) P 


(A3) 


bf^ 


Po' ( ikrBxMQi - Bijtnoj - Bjkkj^^ + BiXpgrmor 


, dniQi . , dboi 

+XpgrBjkj —J + Uijboj + Ujkkj — 


(A4) 
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Po 


ikrBrtnf^ — BijinY^ - Bj^kj 


dm 


,( 0 ) 


+XpgrBjkj- 


dm 


( 0 )\ 


dkr 


+ Uijb^f + Ujkkj 


dkk 

db?^ 


dkk 


+ BiXpgrmf^ 

(A5) 


AXp-uV \2.u?f] (v2 + a/) + f (3w2v + v^) + o?v (v^ + 


15 ( 77 ^ + uP-)^ (P + 


(B 6 ) 


Here mo,, feo, etc. refer to the Fourier space variables for 
simplicity of notation. As a first step in the calculation, 
Eqs. (A2) and (A4) are inserted into Eqs. (A3) and (A5) and 
expanded, neglecting those terms that contain UijUrs, Uij^r, 
QiQj, UijXp, ^Xp’ X% BiBj, BiBij and BtjBrs as higher order 
in this perturbation expansion. 


Appendix B: List of all transport coefficients 



Axj,cP[aP-3P) 
\5p{P + uP)^ 


(B7) 


( 4 «))^ ^ ^ [2aPf,v{u? - 3~P) + iaPr,-^ [P + ca^) 

+2pv(^aP - 3v^) - So)'^ 

X (rj^ + aP^ (P + , (B 8 ) 


In this Appendix we list all transport coefficients 
yS®, ... in the form of integrals over the isotropic 
velocity and magnetic correlation functions, W„ (R, k, to ), 
Hu (k, to) ,Wh (bmR, k, to ), Ht (k, to). This parallels Appendix 
B in RS06 and there is some overlap; however, for complete¬ 
ness we list all coefficients. 

Analogous to the relations in Sec. IV for the Cartesian case 
and RS06, we list here the coefficient of Ank^Wu^i, or AnpHu,h 
in the integrand of each transport coefficient; that is 
or in 


(“hL = j 

f dkdtoPPj^ 

(k, to) Hu,h (k, to ), 


f dkdtoPP'^ 

(k, to) Wu,h (k, (P ), 


An / dkdo) (k, to) Wu,b (k, to). (B1) 


We use the notation 77 = k^t], v = k^v, and V In a = Xag (e-g-, 
Vlnp -H VlnS ^Xpu g) 



XhOiP - lAuPP + 5v 4) 
30p (P + aP)^ 


(B9) 


[ 4 ^^ (1 Iw^V-H -H 477(1 Iw^V-H 

(1 Iw'^v -h Sto^P) + p (l2to^P -P + Uto*) 

- Ap (Sto'^P + 3aPP + 2w®) -H - 5w*j 
X [P + aP)~^ [P + to^Y ^, (BIO) 


ip* ~ + 1 loP^ 

-Apv(P +to^Y h-4^^(-11w‘*v^ -i-5w^v‘* h-Sw®) 

- 8w®v^ H-19w‘*v‘* H-21w*j H- (^P + to^^ p“*, 

(BID 


1. Nonhelical a coefficients 


D’) 

P\ 

(r“D- 

{y'P 


^ Xud 

“ 6 ( 77 ^ -H ’ 


XP 

6 p {P + to^) ’ 


XpuiP 

3 ip + to^) (P + to^) ’ 

^ Xh(^^ - P) 

6p (P + top 


(B2) 


^ [- 477 ^ {3to^v + 5P) - Afi{3to^v + 

- 1 - 77 "^ (44w^v^ -H 13v"^ -H 31w"^) - ^p (3w"^v' -H 

- 2^p (^Sto'^P + 3aPP + 2w®) 

-h5 [%uPP + 3u^P + 

y.iP + aPy^ (v^ H- to^Y^ , (B12) 


(B3) 


(B4) 


(B5) 


-H p (-12w^v^ + iP - 3<y'^) 

- 4?7^ (l7w‘*v^ - 5to^P + lAtP^ 

-56 u^P + 13w'^v‘* - 53w*j (rj^ + to^^ (v^ + ^ p 


-1 


(B13) 



^ [12^2-2-2 ^ ^2^ ^ I2r,^0? - y2) 

+ Au^fjv (\^ + lu?^ + ^if (5w"^v - <y^v^) + Sio'^v'^ - 5w* 

- 77 “* ( 20 w^v^ + 9v"^ + 1 ^ , 

(B14) 

^ [-4w^77v(6w^y^ + y”* + Sw^^) 

- 77"^ (l2w^y^ - 5y"^ + w”*) 

+ 4^^y (6<y^y^ + 5y"^ + 

+ 4/7^ (-3w"^y^ + 3<y^y"^ + 2w®) 

+ 7<y"^y"^ + 9w*j (rj^ + 0 /^ (y^ + cj^'j ^p^*, (B15) 

[_8c^07jv - 1607477'^ - 867^^y 

- 77 "^ (Sw^y^ + y"* + 7(77"^) - 4^^ (7(77"^y^ + 3a7^y‘* + 4a7®) 

+ 12a7®y2 + Sw'^y'* + 7<77*j (77^ + ^ (y^ + , 

(B16) 

(T<"'>), = ^[4iiM-3-V + 2a,=y+3<.‘) 

-Sw^fjv (v^ + +!)“'(—12^3?^ + sy + 

+ 5w''y +11 (y+w-}"'. (B17) 

(r™). = K K - '7+* 7 = (5-.^v+6y) 

+8^ (3(77®y + 4a7^y^) + 16^^ (4a7^y + Sw^y^) 

+ 477 ^ (l3w^y2 + 3u/v‘^ + lOw®) 

+ Scj'^ (4a7^y^ + y"* + 307^^)] 

X (jf + Ci/^ + 67^) , (B18) 

(r<">),-^[-16.7V(3a7V + v^+2a7^) 

+ 77“* (l2(J^ + 19v‘^ - 23a>^) 

- 8 ? 7 ^ (3a7^y + 4a7^y^ + y^) + ff (52ti/v^ + 56li/v^ - 36c<7®) 

+ 40a7®y^ + 37(77"^y‘* - 1367* j (77^ + 07^) (y^ + 67^) p"', 

(B19) 

2. p coefficients 

^ 3 {ff- + 67^) 
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(77"'),, = 0, 

(B21) 

^ 3 {f + 67^) {f + 67^) ’ 

(B22) 

r ,2 . ,2 

M ^ "2’ 

* 6p {f + 67^) 

(B23) 

/jdy)'! - "" 

^ 12(772 + 672f’ 

(B24) 

r,2 , ,2 

(^W) 2’ 

* 12p(y2 + co^') 

(B25) 

2672(11772 - 5 672) 

^ 15 (772 + 672)2 (y2+072)’ 

(B26) 

/ (Q)\ 9y^ - 48672y2 + Jo)'^ 

* 15p(y2 + o72)* 

(B27) 


(«""). = 

774 (2367^ - y2) + I2ff (67^ - 672y2) + 567^ (y2 + of) 

30 ( 77 ^ + 67^)* (y^ + 67^) 

(B28) 

= 

3f [-Uaff + y4 + 367^) - 20674y2 + iSoj^y-i + Uof 
30p {ff + 67^) (y2 + of)^ 

(B29) 

iff^^)^ - ^ |2?7^y + co^) + 1667^?7*y* + 567“* (y^ + co^) 

-i-fj (667‘*y^ - 267®y) - 77 "^ (l067^y^ + 3 y‘* + 767"^) 

- 2^^ (867"^y^ + 367^y"^ + 567®)j + 67^) + 67^) , 

(B30) 

(t**"'). = 15 ("’+“7 - (“'' - 3“'^) 

-Aaffiv (367^y2 + v‘^ + 267“*) - 667®y^ - 67"^y"^ + 367* 

+ ff [-6off + y"^ + 67‘*)j {ff + 67^) {f + of) ^p^', 

(B31) 


(B20) 
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^ ^ [2^'y(5^ + a?) + \ba?r,'^v^ +r,{(^J^v^ - 2u^v) 
+ff (lOwV + 3v^ + Ico'^) + 2f + 3 loH‘^ + 5w®) 


-5u^{v^ + u?)\f + a?)'\y'^+u^\\ 


(B32) 


(;,(«))^ ^ 1 [_4^3-3 (-2 ^ ^ 4^2--^3^2-2 ^ -4 ^ 

+if (-6(y^v^ - 3v^ + Sw"^) + 41;^ [-1(J^\^ - 4w^v"^ + oj^J 

(22w^v^ + Ov"* + w‘*)j + w^) (v^ + p"\ 

(B33) 


3. Helical a coefflcients 



2fl 

3 {ff- + aP-) ’ 


(B34) 


( (0)1 _ 

3p(v2 + w2)’ 


(B35) 


= 0, (B36) 

= 0, (B37) 



?7^ (y^ + 3w^) - cdPv'^ + co^* 
6 ( 77 ^ + (y^ + to^) 


(B38) 


^ f (a;^-y^)-62^(3^ + a7^) 
^ 6 p (i/^ + a>^) (y 2 + 


(B39) 


= _1 [ 3^4 ^^4 _ -4) ^ 4-5 ( 5 ^ 2 - _ 3-3) 

- %u^ifv{^ - 7a/^ + 403“*^y(y^ + 9a/^ 

+ Aif (1 Iw'^y^ + bixP'v' + 5<y®) 

- 5o/ (4w^y^ + v‘^ + 3w"^)j {rf' + ^ + o/'j , 

(B40) 


-Af]^ (3w"^y + 2w^y^ - y^) + 477 ^ (^-lla)‘^\P + 2w^y"^ + 3w®) 

-4^(llw®y+ 18wV +7wV) 

+07“* (-20w^y^ + y"^ + 1 Iw"^)] + w^) (v^ + • 

(B41) 

All of the listed kinematic transport coefficients agree with 
those given in RS06, with one exception. This is the 
coefficient, which contains a factor 1 /30, rather than 1/60. 


Appendix C: The sign of (rjyx)l 

In this appendix we give argue that the sign of (riyx)^ is al¬ 
ways negative, given reasonable assumptions about the form 
of Wb (k, at). We have not been able to find a general proof 
that this is the case due to the complexity of the expression 
Eq. (33), but instead analyze the cases Pm = 1, Pm «: 1, and 
Pm » 1 separately. In addition, plotting (r/yx)^ for Gaussian 
Wb [Eq. (30)] across a range of Pm (e.g.. Pig. 1) leads us to 
the same conclusion for this specific Wb- [Note that (77^x)^ 
depends nontrivially only on Pm and q when written in the di¬ 
mensionless variables given in Eq. (31), meaning it is straight¬ 
forward to observe positivity by plotting {rjyxfb against q over 
a range of Pm.] 


1. Pm = 1 


Inserting v - q into Eq. (33) leads to 


(.Tlyxfb = / doDdkk^Wb {k, u) 


+ u'^) 


An integration by parts in u yields 


i'Uyx)b - 


An 

15 


dojdk 


1 


■ tan 


15 ( 77 ^ H- uP)^ 


07 \ dWb 
77 / du 


5rf -H 307^ dWb 


07- 


{rf- H- 07^) du 


(Cl) 


(C2) 


Under the reasonable assumptions that udWjdu < 0 and 
tan'* (07) dWjdu < 0, each term in the integral must be neg¬ 
ative. (Note that the tan'* (07) 7 /W/( 7 o 7 < 0 condition, al¬ 
though it may appear less familiar, is just as restrictive as 
udWfdu < 0, given the odd nature of the tan'* function.) 


2. Pm <K 1 

Inserting 77 = y/Pm into Eq. (33), we carry out a series ex¬ 
pansion about Pm * = 00 of the resulting expression. The 
reason for this expansion (rather than the more obvious ex¬ 
pansion about Pm = 0) is that we wish to explore that low Pm 
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limit with large i] rather than that with v —> 0, since SOCA term is obviously negative due to the positive definiteness of 

looses applicability as v,ri 0. The series expansion to first the integrand, 

order in 1 /Pm“* is 




J dcjdkWbk^ 

4v2 




• lay 


(v^ + OJ^)^ 
1 


15 (v2 + w2) Pm 


+ ... 


(C3) 

(C4) 


The first term is independent of Pm, persisting as ^ 0, and 
the existence of this is not surprising given the fact that the 
dynamo can arise from the B Vb + b- VB term in the induction 
equation. This term can be shown to be negative using the 
same integration by parts method used to obtain Eq. (C2), with 
the requirement udWjdu < 0. The Pm dependent second 


3. Pm » 1 


Inserting v = Pm 77 into Eq. (33), and carrying out a series 
expansion about Pm = 00 (see previous section), one obtains 


idyx)h ~ 1C 


\6n 

15 


dcjdk WhK 


' 1 K-f) 

^ Pm (jj 2 _|. ^ 2^2 



(C5) 


As expected, there is no v = 0 contribution to the transport. 
Again using integration by parts, it is easy to prove negativity 
of the integral provided codWIdco < 0. 



